Ehrenfest-time dependence of weak localization in open quantum dots 
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Scmiclassical theory predicts that the weak localization correction to the conductance of a bal- 
listic chaotic cavity is suppressed if the Ehrenfest time exceeds the dwell time in the cavity [I. L. 
Aleiner and A. I. Larkin, Phys. Rev. B 54, 14424 (1996)]. We report numerical simulations of weak 
localization in the open quantum kicked rotator that confirm this prediction. Our results disagree 
with the 'effective random matrix theory' of transport through ballistic chaotic cavities. 
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The wave nature of electrons is the cause of striking 
effects that are absent in classical mechanics. The best 
known manifestations of quantum mechanics on electri- 
cal transport are weak localization (WL), universal con- 
ductance fluctuations (UCF), and shot noise In this 
letter we consider such quantum mechanical effects in a 
ballistic cavity with chaotic classical dynamics. A phe- 
nomenological theory of shot noise, weak localization, 
and conductance fluctuations in a ballistic chaotic cavity 
is given by random matrix theory 0] . The predictions of 
random matrix theory have been confirmed using semi- 
classical methods [3, y, |3, |(j . 

Since they are intrinsically quantum effects, it is legit- 
imate to ask what happens to shot noise, weak localiza- 
tion, and conductance fluctuations if the classical limit 
is taken. The classical limit corresponds to the case that 
the electron wavelength is much smaller than the system 
size. The relevant time scale for this question is the so- 
called "Ehrenfest time" 0, 0} , the minimal time at which 
the wave nature of electrons become apparent. This is 
the time at which a minimal wavepacket is stretched to 
the system size L (or the size of the lead opening) . It is 
given by 



te = A In kL, 



(1) 



where A is the Lyapunov exponent characterizing the 
chaotic motion in the cavity, and k the electron wavenum- 
ber. 

The effects of the Ehrenfest time are most prominent 
if Tg is larger than the mean dwell time in the cavity td- 
In this regime most of the transport is classical and shot 
noise is suppressed |a 111 LL0| • The effect of te on WL was 
first addressed by Aleiner and Larkin Q , who predicted a 
suppression proportional to exp(— 2te/td). Exponential 
suppression of WL, but with a different exponent, was 
also found by Adagideli . Suppression of WL was ob- 
served experimentally by Yevtushenko et al. fl^ ] . There 
is no semiclassical theory of the Ehrenfest-time depen- 
dence of UCF. 

Numerical simulation of systems with large Ehrenfest 
times is difficult due to the need for (exponentially) high 
values of kL. To reduce numerical costs Jacquod et al. 
proposed to replace the two-dimensional cavity by a one- 




FIG. 1: Schematic picture of a ballistic chaotic cavity of 
size L. The cavity is attached to two leads, labeled "L" and 
"R", and with width d^ and drt, respectively. The contours 
Cl and Cr drawn in the leads are used for the calculation of 
conductance. A schematic drawing of relevant trajectories for 
weak localization is depicted inside the cavity. (The physical 
trajectories scatter from the cavity walls many times.) 



dimensional quantum map |13| |. The map is 'opened', 
so that simulation of transport properties is possible. An 
example of such a map, the open kicked rotator, was used 
to obtain numerical results for shot noise WL (lif . 
and UCF HHHGl. The simulation results for shot 
noise were in agreement with semiclassical predictions Q . 
In contrast, the results for WL and UCF are remarkably 
different: no dependence on te was found. To explain 
their results, the authors of Refs. [H, EE [H [3 pro- 
posed an 'effective random matrix theory': Noting that 
quantum diffraction contributes only if the classical dwell 
time is larger than te, they proposed to apply random 
matrix theory only to the part of phase space with dwell 
times longer than te 01- Since UCF and WL are not 
dependent on the size of the phase space involved, the 
'effective random matrix' theory is able to describe the 
numerical results for WL and UCF 20]. 

Apart from the prediction of the size of the weak local- 
ization correction, the 'effective random matrix theory' 
and the semiclassical theory have different predictions of 
the minimal time required for weak localization to oc- 
cur. In the semiclassical theory, quantum interference 
requires a minimal wavepacket to be split and reunited, 
which takes a minimal time 2te- A schematic diagram 
drawing relevant semiclassical trajectories for WL is de- 
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picted in Fig. ^ This is different in the 'effective random 
matrix theory, where quantum interference is fully estab- 
lished already after a time te- 

It is the purpose of this letter to resolve the discrep- 
ancy between the reported numerical simulations and the 
semiclassical theory for weak localization. In the first 
part of this letter, we review the original semiclassical 
theory Q of the Ehrenfest time dependence of WL and 
amend it in order to correctly implement classical correla- 
tions in an open cavity [21| . In the second half we report 
time-resolved numerical simulations of weak localization 
in the open quantum kicked rotator. The time-resolved 
simulations allow us to determine the time at which WL 
first appears, as well as the effect of a finite Ehrenfest 
time on the conductance. We find that the time-resolved 
simulations for WL and the corrected semiclassical the- 
ory are consistent. 

Semiclassical theory. The semiclassical theory of Refill 
considers a two-dimensional chaotic cavity which is cou- 
pled to two electron reservoirs via contacts of widths 
(Jl,R) see Fig. One has c?l,r <C L, so that electron 
dynamics inside the cavity is ergodic. The system is stud- 
ied in the semiclassical limit, where the number of open 
channels is large, c?l,r& 3> 1. This separation of length 
scales is equivalent to the following separation of time 
scales (vk)^ 1 <C L/v ~ 1/A <C td, where v is the elec- 
tron velocity, and the dwell time is td = nA./(di, +dn)v, 
A being the area of the cavity. It will be useful to define 
the probability that an electron at a random point will 
leave the cavity through contact i, Pi = di/(dji + di). 

The central objects in the semiclassical theory of 
Aleiner and Larkin are the non-oscillating parts of the 
product of advanced and retarded Green functions 
The "diffuson" V is composed of combinations of orbits 
with themselves while the "cooperon" C is composed of 
orbits and their time reversed counterparts. Quantum 
effects are introduced by including a weak additional ran- 
dom quantum potential. This random potential leads to 
two different effects: phase space diffusion, and coupling 
between paths which are close in phase space. Both are 
expected due to quantum uncertainty. After averaging 
over the random potential, the equation of motion for 
the leading order diffuson, T>°, is 



V°(u; 1,2) =6(1,2) 



(2) 



where j = (0j,Rj) = 1,2 denotes the phase space co- 
ordinates, limited to the energy shell: R is the elec- 
tron position while the angle 4> is the direction of its 
velocity v. The strength of the phase space diffusion is 
7 q ~ X 2 /kv. The operator L\ is the Liouville operator. 
For a hard wall cavity, L\ = v% ■ d/dTLi, with appro- 
priate boundary conditions at the walls. The symbol 
(5(1,2) = 27T(5(Ri — R,2)S((f>i — 4>2) denotes a delta func- 
tion on the energy shell. The leading order cooperon, C°, 
also satisfies Eq. ©. 



The WL correction to the diffuson arises from combi- 
nations of a trajectory with small angle self-intersection 
and a partner without such self-intersection 0, , see 
Fig. n The 'interference region' is composed of close 
phase space points (with distance ~ fc -1 / 2 ) where quan- 
tum diffraction is important. This leads to a cjuantum 
correction for the diffuson V = T>° + AT> with 



AP(1, 2) = P°(l, 2)^%^ + 2) 



2nTiv 
d3P°(l,3)X>°(3,2) 



2-kUv 

d 2 1 C°(3,3) 
2irhi> ' 



dcj> 2 



(3) 



where j = (cj)j + tt, Rj) denotes the time reversal of phase 
space point j and v = m/2Trh 2 is the density of states per 
unit area. Since we consider DC transport, the frequency 
lo has been set to zero. 

The total (DC) conductance of the dot can be ex- 
pressed in terms of the diffuson T> 0, 0] , 

G = 2vtivv 2 dh -p. cos(0 R - 0!) 



dl 2 



0L+7T/2 



C L J<Pl-7v/2 



2lT 



cos(0 L -<fe)£>(l,2). (4) 



Here G is measured in units of the conductance quantum 
2e 2 /h, and Cl and Cr denote cross section contours of 
the contacts. The angle 0l is the direction of the inward- 
pointing normal to Cl , while 4>R is the outward-pointing 
normal to Cr. 

The WL correction to the conductance 6G is obtained 
by substituting Eq. © for AV into Eq. gj. The first 
two terms in Eq. (jSJ) contain paths which leave the dot 
at the point of entry, and hence do not contribute to the 
transmission. The phase space point 3 in the third term 
in Eq. can be viewed as the center of the 'interference 
region' depicted in Fig. ^ As the two diffusons and the 
cooperon head into opposite directions in configuration 
space, one can treat them as statistically uncorrelated 
and estimate them independently 0|. 

Both the cooperon and the product of diffusons in Eq. 
exhibit non-trivial correlations 0], since the incom- 
ing path is strongly correlated with the outgoing one, cf. 
Fig. ^ An intuitive way to estimate the cooperon (or 
the product of diffusons) was presented by Vavilov and 
Larkin . They demonstrated that the effect of phase 
space diffusion is equivalent to an average over phase 
space regions of size ~ k~ x / 2 near 3 and 3. Following 
Ref. 12.1 we consider phase space points 3' and 3" within 
a distance of order fc -1 / 2 from 3. The phase space dis- 
tance between classical orbits starting at the points 3', 
3" will diverge exponentially due to the chaotic dynam- 
ics and will be macroscopic after a time te/2. Once the 
phase space distance between the trajectories is large, 
they can be considered uncorrelated. 
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For the cooperon, the classical correlations mean that 
it takes a minimal time to form a closed loop. Let us 
denote by tj the time it takes for a classical trajectory 
starting at phase space point j to leave the cavity. If 
£3 < Te/2 all the points near 3 leave the cavity before it 
is possible to close a loop. In contrast, when £3 > te/2, 
the cooperon is estimated by propagating phase space 
points close to 3 forward in time, and points close to 3 
backward in time, for a time of te/2. After this propaga- 
tion the points which contribute to the local average are 
uncorrelated. Then the cooperon can be estimated using 
its ensemble average m/A. The resulting cooperon in 
the open cavity is given by 

C°(3, 3) = <C°(3", 3"')) * ^9(t- 3 - r E /2), (5) 

where 8(x) = 1 if x > and otherwise. 

The product of diffusons is estimated similarly. The 
phase space points 1 and 2 in Eq. (|3J) refer to two different 
contacts. If t% < te/2, trajectories starting close to 3 are 
still correlated when they leave the cavity. Therefore, 
they cannot arrive at points on different contacts, leading 
to a vanishing contribution to SG. When t$ > te/2 the 
lead integrals will result in a factor of PrPl, as in Refs. 
0, El Note that the WL paths are composed of four 'legs' 
with minimal time of te/2, so that the minimal time of 
the whole path is 2te- 

The resulting WL correction is 

SG = -PrPl^J J d30(t 3 - r E /2)L 3 9(t 3 - r E /2). (6) 

The Liouville operator in Eq. Jfjjl measures the rate of 
flow of probability density out of the integration range 
of the phase space variable 3 The boundary of the 
integration range is composed of the lead phase space 
points propagated backward for a time te/2. This leaves 
only a fraction exp(— of the size of the boundary 
at the lead. In addition, only points which also have 
t 3 > Tg/2 will have a non- vanishing cooperon, leading to 
another factor of exp(— ^ E -). Hence, we find 

6G=-P K P L e-™/™. (7) 

The result @ shows that, even after accounting for clas- 
sical correlations that were not accounted for in Ref. 0, 
the WL correction to transmission is suppressed expo- 
nentially. However, the exponent is different from that 
of Ref. 3| where it is reported that SG oc exp(— 2te/td). 

Numerical simulations. We have compared the predic- 
tions of the semiclassical theory to numerical simulations 
of the open quantum kicked rotator. The time evolution 
of the kicked rotator is discrete, and given in terms of 
the Floquet operator ip(t + l) = J-tp(t). The state of the 
system ip is a finite vector of dimension M. The Floquet 
operator is an M x M matrix and reads 



where n is a numerical parameter that determines the 
classical dynamics of the map and < 9 < 2tt/M a 
parameter that sets the precise quantization condition. 
The region k > 7.5 is associated with classically chaotic 
dynamics. The size of the matrix M is taken to be even. 

To open the system, two consecutive sets of N elements 
of ip are chosen as the contacts 01 • A scattering matrix 
is then constructed by the rule 

00 

S(e) =J2e iet S(t); S(t) = P [TQ] t ^ 1 TP T , (9) 
t=i 

where e is the quasienergy, P is a 2N x M matrix pro- 
jecting on the contact sites, and Q = 1 — P T P. The 
semiclassical limit corresponds to M — > 00 while keeping 
td = M/2N fixed. For large values of k the Lyapunov 
exponent A can be approximated A ~ In (k/2), yielding 



TE 



In TV 
ln( K /2) : 



(10) 



up to an A-independent constant. Once the scattering 
matrix is known, the conductance G follows from the 



Landauer formula, G(e) = J2 m =i J2n=N+i \Smn(e)\ 2 - 

We study transport properties as a function of time, 
simply by truncating the sum in Eq. @ after a time 
t = to. We calculate the ensemble averaged conductance 
by averaging over e, 9, lead positions, and a small range 
of k. An advantage of the time-resolved approach is the 
ability to evaluate the quasienergy average analytically. 
In order to achieve sufficient accuracy, we take between 
20 000 and 80 000 realizations for 9, the lead positions, 
and k. 

Ideally, for two contacts with N channels, the ensemble 
average of the classical conductance is N/2 and WL can 
be calculated by taking the difference between the quan- 
tum and classical conductances. In practice, however, 
the ensemble average of the classical conductance may 
differ from N/2 and this simple procedure fails. Alter- 
natively, one compares conductances with and without a 
time-reversal symmetry breaking perturbation. However, 
such a perturbation may affect the classical dynamics, 
complicating the interpretation of the results. Instead, 
we choose to subtract the classical conductance by com- 
paring conductances at N and 27V channels, for the same 
realization of the classical parameters. The classical con- 
ductance is proportional to the channel number, whereas 
the WL correction is independent of it. 

Simulation results for the WL correction to the conduc- 
tance for td = 5 and k ~ 10 are shown in Fig.JSJas a func- 
tion of the cut-off time to- We also performed simulations 
for td = 5 and n w 20, td = 10 and K « 10, and td = 10 
and K s» 20. In all cases, the derivative dr^/dlnN, as 
extracted from the low-to part of the curves, is consistent 
with an onset of the WL correction after a time 2te, in 
agreement with the predictions of the semiclassical the- 
ory. A summary of our results for SG is shown in the 



^2N 



F„ 



(8) 



x e 



-i(KM/47r)(cos(27rm/M+0)+cos(27rn/Af+0)) 
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FIG. 2: Time-resolved weak localization correction SG as a 
function of the cut-off time to, for k ~ 10 and to = 5, and 
for four values of N. Inset: summary of results for the weak 
localization correction. Data points are obtained from the 
limit to to of simulations for to — 5 and n ~ 10 (+), 
to ~ 5 and k ~ 20 (x), to = 10 and k ~ 10 (*), and to = 10 
and k ~ 20 (o). The solid lines are fits oc exp(— Te/ttj), with 
dTE/dlnN taken from the 1ow-£q part of the simulation data. 



inset of Fig. [3 together with a fit oc exp(— te/td), with 
drs/dlnN obtained from the low-io part of the data. 
The simulated SG has a systematic dependence on te, 
consistent with the semiclassical theory. Previous nu- 
merical calculations of SG did not show a dependence on 
te probably due to insufficient accuracy. The fact 
that the over-all magnitude of WL is smaller than in Eq. 
(J7J is well known |l5j and caused by a lack of ergodicity 
for the dwell times we consider. Note that our results for 
both the onset of WL and its magnitude disagree with 
the 'effective random matrix theory'. We have further 
verified our results by calculating SG for the three-kick 
model, a version of the open quantum kicked rotator in 
which time-reversal symmetry can be broken jl5j . We 
found SG = when the symmetry is broken. 

We believe that our simulation results, together with 
the corrected semiclassical theory, resolves the contro- 
versy around the Ehrenfest-time dependence of the weak 
localization correction to the conductance. Our results 
quantitatively confirm the corrected semiclassical theory 
and disagree with the 'effective random matrix theory'. 
The results presented here are limited to weak localiza- 
tion and do not carry over to conductance fluctuations. 
We have also carried out numerical simulations of the 
variance of the conductance for the open quantum kicked 
rotator |2^|, and find behavior consistent with that re- 
ported previously in the literature 0, 0, 0: no de- 
pendence of conductance fluctuations on the Ehrenfest 
time te and an onset time te, not 2te- Currently, there 
is no semiclassical theory these results can be compared 
with. 
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